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On q-analogues of Riemann's zeta 

By Ivan Cherednik* 



The aim of the paper is to define q-deformations of the Riemann zeta 
function and to extend them to the whole complex plane. The construction is 
directly related to the recent difference generalization of the Harish-Chandra 
theory of zonal spherical functions [C1,C2,C3]. The Macdonald truncated 
theta function [Ml] replaces x s . The analytic continuation is based on the 
shift operator technique and a more traditional approach close to Riemann's 
first proof of the functional equation (see [E,R]). 

We introduce g-deformations using the real and imaginary integration. 
The construction depends on the particular choice of the path and the inte- 
grand much more than in the classical theory. We come to several different 
functions with different properties. All have meromorphic continuations to 
the whole complex plane. The basic ones approach £ (up to a T-factor) as 
q — ► 1 for all s thanks to the Stirling-Moak formula for T q [Mo] . There are also 
equally interesting transitional q-zeta functions converging to £ for 5Rs > 1 and 
to proper combinations of gamma functions for 3fts < 1. They approximate £ 
very well in the critical strip when q is close enough but not too close to 1. 
Then the tendency T q — ► V suppresses the Stirling formula and eventually they 
go to their T-limits. 

Our q- integrals do not satisfy the celebrated relation between and 
C(l — s) in the g-setting as well as the Euler product formula. However certain 
analytic properties of the g-zeta functions are better than the classical ones. 
For instance, the imaginary one is periodic in the imaginary direction and 
stable as Sfts — * ±oo. The number of its zeros can be calculated exactly upon 
the symmetrization s <-> 1 — s. Its leading term coincides with that from the 
Riemann estimate of the number of zeros of £(s) in the critical strip. The 
g-zeros do not belong to the critical line anymore but their distribution is far 
from random. 

Let us try to summarize the key results of the paper in more detail. 
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1) We construct the meromorphic continuation of the imaginary g-zeta 
functions (defined in terms of the imaginary (/-integration) to all s by means 
of the shift operator method, mainly Theorem 4.2. The existence of such a 
continuation is not too surprising because we can also apply Cauchy's theo- 
rem almost following Riemann's first proof of the functional equation. The 
functions with poles dense everywhere appear in this approach. 

2) What is surprising is that the convergence of the meromorphic continua- 
tion to the classical zeta (up to a factor) holds for the g-plus-zeta corresponding 
to the integration with the kernel (q x + 1) _1 , but fails for the one defined for 
(q x - The latter tends to T(s)((s) for 9fo > 1 only. In the plus case, the 
limit for negative Sfts is calculated via the shift operator. We were not able to 
empoy other methods. 

3) Numerical experiments show that the g-plus-zeta does not have zeros 
for 9fts > 1/2 when < q < 0.96. Actually it is very surprising because the 
imaginary q-zeta functions are periodic in the imaginary direction and the 
convergence to the classical zeta (multiplied by a proper factor) is for relatively 
small 9s. Presumably all zeros of the g-plus-zeta come from the classical zeros 
for > and a > 1. 

4) The imaginary g-integral can be exactly calculated for the kernel q~ x . 
It is an important part of the paper (Section 2) related to Ramanujan's q- 
integrals and closely connected with [CI]. The resulting function has no zeros 
for ?R.s > 1/2. However it explains the above numerical phenomenon for small 
q only. Analogous exact formulas are obtained for the Jackson summation 
instead of the imaginary integration. 

5) The most promising result of the paper may be the definition of the 
sharp g-zeta functions, which have no classical counterparts. They are Jackson 
sums in terms of the plus-minus kernels discussed above, convergent for all s 
except for the poles. Their limits q — > 1 are similar to those for the imaginary 
integration. The analysis of their zeros leads to some interesting properties of 
Riemann's zeta. 

6) Still the main applications are expected upon the (anti)symmetrization 
of the integral g-zeta functions with respect to s <-> 1 — s (Section 6). In the 
imaginary case, the number of the corresponding zeros matches well the famous 
Riemann estimate. The antisymmetrization of the sharp g-zeta may satisfy the 
Riemann hypothesis. Conjecturally, all its zeros in a natural horizontal strip 
are imaginary. 

Actually our methods are of more general nature. They can be applied 

2 

to study difference Fourier transforms of functions of q x -decay for a quadratic 
form x 2 associated to a root system in R™ (or in the imaginary variant). 
Without going into detail, we mention that such class of functions is Fourier- 
invariant in the g-theory in contrast to the Harish-Chandra transform. The 
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q- Fourier transforms "at zero" of (q~ x ±1) _1 , simplest deformations of the q x , 
are natural multi-dimensional counterparts of the (/-integrals studied below. 

Note that q-zet& functions are mainly considered in the paper as a good 
starting point for the general analytic theory of the (/-Fourier transform. We 
use them to demonstrate the shift operator method and to show some typical 
problems. Thanks to such a concrete choice (and the one-dimensional setup) 
we are able to illustrate all claims numerically. There are of course more 
specific results and applications to the classical zeta, which hopefully will be 
continued. 

The first section contains a certain motivation and the main constructions 
(for the standard minus- kernel) . Section 2 is devoted to the exact formulas for 
the integrals of the truncated theta with and without the Gaussian. In Section 
3, we introduce the (/-zeta functions systematically using the real and imaginary 
integrations and the analytic continuation. Section 4 contains the approach 
via the shift operator. The next one is a discussion of the (/-zeros. The last 
section is an attempt to establish connections with the Riemann hypothesis. 

The author thanks P. di Francesco and E. Opdam for useful discussions. 
I acknowledge my special indebtedness to David Kazhdan, who greatly stimu- 
lated this work. 



1. Main constructions 

Switching from integers to their (/-analogues is a common way of defining 
g-counterparts. For instance, one can deform £(s) = J2T=i n ~ S as follows [UN]: 

(1 „ <f( . ) = E ^ faWf tfi=i.. 

n=l 1 q 

Here we assume that < q < 1 and > 0. The numerator ensures the 
convergence. It becomes a constant when Q q is rewritten in terms of 
This function approaches ( as q — ► 1 for 5Rs > 1. We note that similar q- 
deformations were considered before in connection with Carlitz's Bernoulli 
numbers and p-adic L-functions (see [Sa]). Ueno and Nishizawa constructed a 
meromorphic continuation of ( q to all s with simple poles in { — Z + + 2iraiZ} 
for a = — (log((/)) _1 , Z + = {0 < n € Z}. The procedure is a certain substitute 
of the functional equation which cannot be saved. As a rule relations between 
special functions which involve the multiplicative property of x s have no direct 
^-counterparts. 

A modification of Cq tends to £ for all s. The latter property is impor- 
tant but not sufficient. Once the functional equation is sacrificed, something 
else should be gained. At least, certain analytic advantages can be expected. 
Generally speaking, ^-functions (and their representations as sums, products 
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or integrals) are single-valued and have better analytic properties than their 
classical counterparts. It does not seem to happen with ( q (s). Actually its 
definition is a transition from the differential theory to the difference one be- 
cause the complex powers are still involved. The latter are redefined via T q in 
a systematic g-theory. 

Indeed, substituting q sn — ► q sn / 2 in (1.1), the resulting zeta (up to a 
simple factor) is the classical ((s) where x s is replaced by sinh s (x/(2a)). By the 
way, such numerators are more relevant to associate Cq( s ) with the quantum 
SL2. This passage is well-known in the theory of special functions. Instead 
of the Bessel functions, which are pairwise orthogonal with respect to the 
measure x s dx, one gets the hypergeometric and spherical functions. The multi- 
dimensional theory of sinh(x) s is the cornerstone of the harmonic analysis on 
symmetric spaces (see e.g. [He,HO]). It is of obvious importance to incorporate 
the hypergeometric function into the theory of zeta. The paper [UN] is a step 
in this direction. However the basic (difference) hypergeometric function seems 
more appropriate. 

A straightforward definition. Let us switch from [x] 2k to the following 
combination of g-shifted factorials (see e.g. [A,M1]): 

( L2 ) S k (x;q) = [I ^ _ _ qj+k -2xy 

Trying to improve the ze q due to Ueno - Nishizawa, we put 

c 

Q(s) c = ^L S (M?), Q(s) = gOOoo for 

I 1 - 3 ) n=l 

q = exp(— 1/a), R 3 a > 0, u G C, c = 1, 2, . . . . 

From now on, q and a will be always connected as in (1.3). Since 5(x) is even 
the sign of can be arbitrary. As we will see, y/a and the "direction" u are 
necessary to ensure the convergence and go back to the classical theory via the 
following Stirling- Moak formula [Mo, UN]. 

Lemma 1.1. a) Setting z k = exp(/clogz) for complex {fixed) z, k and the 
usual branch of 'log, 

lim (a/4) fc <5 fc ( v / oI;g) = {-z) k 
(1.4) a ^°° 

if either z R + or k G Z. 



b) Replacing az by v for any (fixed) v G C 
(1.5) lim(a) 2k 5 k (V^;q) 



2kx , ^..^ _ T(k + 2^)T(k-2^) 



r(2 v / ^)r(-2 v ^) 
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The singularities of the denominator in (1.4) go to infinity together with a 
for z G" R + and disappear in the limit. When k G Z the claim is obvious since 
<5fc is a rational function in this case. The Moak theorem also gives a certain 
estimate of the convergence. The second formula is well-known. 

It is clear why the first formula fails for z G R+. The numerator of (1.2) 
becomes for infinitely many z provided that k G" Z. However the limit cannot 
be zero. Another reason is that we switch to (1.5) when z = v/a for fixed v. 
The a-factor in this formula is a 2k = o(a k ) for $lk < 0. Hence (1.4) diverges 
as a — ► oo for z approaching zero and such k. The latter argument will be 
applied to analyze the limiting behavior of the sharp q-zeta functions. 

Thus the simplest choice u = 1 in (1.3) is impossible, since it destroys the 
termwise convergence to the classical sum for ((s). The lemma gives that 

c 

\im(-Au/a) s Q(s) c = £ n~ s 

n=l 

for any s G C (and fixed c) if u R+. However we must avoid u G — R+ 
because lim^oo Cq( s )c does not exist for such u. Let us examine the imaginary 
unit u = i. 

Theorem 1.2. Let ( q (s) = Q(s), 5Rs > 0. Then ( q (s) exists apart from 
the set S(a) of the poles of all terms in (1.3). If $ts > 1, given a sequence 
a m — ► oo and k U^S(a m ), 

oo 

(1.6) lim (-4i/a m ) s ( q (s) = C(s) = 

m— >oo • 

n=l 



We see the first special feature of the g-zeta functions. The convergence 
is for -Rs > in contrast to the usual inequality Sfts > 1. However (1.6) holds 
for Sfts > 1 only. 

The function ( q (s) is not real for real s. The following function is: 

oo 

(i.7) ^0Mv^ ; v) + ^((V»);?))/2. 

n=l 

Here and further z* is the complex conjugation. The limit of this function as 
a — > oo remains the same. 

The analytic behavior of Cq( s ) 1S totally singular. The set S(a) is dense 
everywhere in C unless 2ira G Q. If 2ira = n/d for coprime integers n,d, then 
S(a) belongs to the union Q of the translations of the diagonal cross C = 
{x 2 G iH} by {(l/d)i, I G Z}. Choosing a m = m 2 a, we have the embeddings 
S(a m ) C S(a) and can apply the theorem at least when s Q. In general, the 
sum for ( q (s) is pointwise convergent in spite of infinitely many poles in any 
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neighborhood since the residues are getting smaller faster than the distances 
to any fixed s S(a). 

It seems that "totally singular" functions are inevitable in the difference 
calculus. However we prefer to avoid them in the basic definitions. The next 
step is to remove as many poles as possible turning from sums to integrals. 
We note that ( q (s) will be directly related to 3+ q (k\d) from (3.25) playing 
important role in the analytic continuation of the imaginary g-zeta functions. 

Real integration. Actually the 5-function and its multi-dimensional 
generalizations due to Macdonald are important because they make selfadjoint 
remarkable difference operators which have many applications (see [M2,C3]). 
Hence it is more logical to use them as measures, which is exactly what we are 
going to do. 

From now on we will switch from s to the "root multiplicity" k = s — 1/2, 
standard in the harmonic analysis on symmetric spaces and in the correspond- 
ing theories of orthogonal polynomials. This choice is convenient in the theory 
of ( as well: the critical line 3?s = 1/2 coincides with the fc-imaginary axis. 

First, we will express the T- function and (F in terms of the Gaussian 
exp(— x 2 ) modifying a little the classical representations: 



Since the first formula has a nice ^-generalization with 5k(x;q) instead of x 2k 
(see [CI] and the next section), we can try to use the second to deform 3(&)- 
Let 



The contour of integration C± for 3 goes from to ±ez and then to oo ± ei 
for fixed e > 0, which is assumed to be sufficiently small. Note that we cannot 
integrate over R + for real k because <5& has infinitely many poles there. In 
contrast to (1.8), there is no problem with the point x = 0. Indeed, 5^ has 
a zero of the second order at which annihilates the pole of the same order 
coming from (exp(x 2 /a) — l) -1 . To be more exact, we need to exclude k = 
in this argument. The functions 3^ are regular apart from 



(1.10) K e ± (a) = {2C e ± - Z+ + {2itai)Z} U {-2C e ± - Z+ + (2vrai)Z}. 



They are obviously 27rai-periodic. However it is of no use, since the connected 
components of C\K±(a) have the same periodicity. In this definition, we loose 
the reality of £ on the real axis. The following functions have this important 



(1.8) 




(1.9) 
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property: 

let 3;(k)+t+(k) _ 1 /•+»-< S k (x;q) dx 



(1.11) w s5M = ;/_ 



oo — e« 



exp(x 2 /a) — 1 ' 



(112) 3{ w g V»)-3fW = ^ 



/oo+« exp(x 2 /a) - 1' 

The second integration is with respect to the path {Ci — C e + }. The corre- 
sponding function is well-defined in the horizontal strip {— 2e < 9fc < +2e} 
and has no singularities for 3ft A; > 0. The natural domain of the first function 
is K re d = f C \ {K £ _ U K%}. 

Theorem 1.3. a) The functions 3^ coincide when 

(1.13) 3& + 27raZ [-2e,+2e]. 

XTie functions 3™, 3g /iaue meromorphic continuations from R + to C. T/ie /irst 
/ias no po/es in t/ie strip 

(1.14) if a (a) d = f {k, -2e a < 9fc < 2e a } for e a = 7 min{V^, vra/2}- 

TTie set of poles of the continuation o/3q in strip is {— Z + /2}. It a/so has 
simple zeros at k = 1, 2, . . . . 

&) Let 3(ife) = C(k + l/2)T(k + 1/2) (as above) and $lk > 1/2. Then 

(1.15) lim (a/4) fc - 1/2 exp(±7rifc)3^(A0 = 3(*0- 

(1.16) lim (a/4) fc_1 / 2 3S(fc) = sin(mfc)3(£0- 
For any fegC, 

(1.17) lim (a/Af-^y^k) = cos(^ k)d(k). 



The calculation of the limits of 3^ is based directly on (1.8) and Lemma 
1.1; (1-16) is a formal corollary of the (1.15). 

The limit of 3q € (k) is calculated using the following well-known integral 
representation for the zeta (see e.g. [E]): 

(1.18) (1/2)6 , ;_/ " = -sin(7rg)3(fc) = cos(7rfc)3(fc). 
Joo-ei iyj~z(e z - 1) 

The path of integration begins at z = —ei + oo, moves to the left down the 
postive real axis till — ei, then circles the origin and returns up the positive real 
axis to ei + oo (for small e). Note the difference from the integral f™7^ above, 
where the orientation is different and the path goes through the origin. The 
left-hand side multiplied by (a/4) fe_1 / 2 is exactly the integral (1.11) after the 
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substitution: z = x 2 /a. The factor i\J —z is {a/2)dz/dx, \fx is the standard 
branch positive on R + . 

In the paper, we prefer to replace (exp(x 2 /a) — l)^ 1 by (exp(x 2 /a) + 
1) _1 , which improves the analytic properties of the g-zeta functions. For in- 
stance, the left-hand side of (1.16) diverges for $lk < 1/2, whereas its "plus- 
counterpart" approaches the zeta (up to a proper factor) for all k. However 
the asymptotic behavior of 3>\{k) for $lk < 1/2 is very interesting: 

(1.19) lim o 2fc_1 3*(A;) = tan(mfc)r(fc) 2 for m < 1/2, 

a— too H 

(1.20) {a/Af- 1 ' 2 ^) = (4a) 1 / 2 - fc tan(^)r(fc) 2 + sin(mfc)3(£0 + o(l) 

for ?R.k = 1/2, where the error o(l) tends to as a — ► oo. In fact the latter 
formula holds for all k. Indeed, the term sin(7r/c)3(fc) is suppressed in (1.20) 
for ^tk > 1/2 and vice versa for K/c < 1/2 because of the difference in the 
a-factors. Moreover, the error is of order tan(7rfc)r(/c) 2 in the strip — 1/2 + e < 
Jifc < 1/2 — e for any small e > 0, when a are big but not too big. Since 
r(/c) 2 quickly approaches as Qk increases, the function (a/A) k ~ l / 2 ^\{k) first 
approximates sin(7rA:)3(fe) with high accuracy and then slowly switches to the 
the leading term (4a) 1 / 2 ~ fc tan(7r/c)r(/c) 2 in the critical strip. Numerically, the 
latter (dominant) regime is very difficult to reach for big $sk, even for k ~ lOi. 

The main advantage of 3>\(k) is the following representation (Cauchy's 
theorem), which has no counterpart in the classical theory: 

*Uk) = a7T f\ (l-e-^ +fc )/ B )(l-e-fa--*)/") 
°<?W 2 11 n - e -Cj+2fc)/a)M _ e -(i+i)/a) * 

j—Q \ I \ ) 

~ X (1 - e -(l+2k-l)/ a \ (1 _ Jl+k)/a\ I 
(121) ^PTT { ,[ ' - 

hti (1 - ^ (/+ ^ 1)/a )(! - e ' /a ) exp^)-!' 

Thanks to the factor (exp((/c + Z) 2 /(4a)) — 1) _1 ~ exp(— Z 2 /(4a)), it converges 
for all k (0 < q < 1) except for the singularities, and very quickly. For instance, 
we can caluclate 3q(k) numerically for a ~ 10000 (i.e. when q ~ 0.9999), and 
Qk ~ 100. For integral k, the range is almost unlimited. The numerical 
simulation of the 3 r q e (k) is approximately 10 • • • 100 times worse. 

In a sense, the series (1.21) is complimentary to Cq( s ) (Theorem 1.2). They 
appear together in the following construction. 

Imaginary integration. The integration in the real and imaginary direc- 
tions are practically interchangable in the classical theory of zeta function (the 
substitution x — ► ix multiplies x s by a constant). In the harmonic analysis, 
the difference is fundamental. The imaginary case (integrating with sin(x) 2fc ) 
is much simpler than the real one (for sinh(x) 2fc as the kernel). The real and 
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imaginary integrations are closer to each other for 5k(x; q). However the latter 
has certain analytic advantages. 
For Wt > 0, we set 

(1.22) 3f(A:) = f (-i) / (ex P (-x 2 /a)-ir 1 5 k (x;q)dx. 

J o 

The complex conjugation is denoted by *. 

Theorem 1.4. a) The function 3 l ™(k) is analytic for $tk > 0, and has a 
2irai-periodic meromorphic continuation 3q W (k) to all k £ C mi/i the following 
set of poles (simple for generic a): 

(1.23) V = {-Z + /2 + iraiZ] U A, /or 



A = {-2v / 2vraiN - Z+ + 2vraiZ} U {-2v / 2vraiN - Z + + 2vraiZ}. 

Tae poZes o/ the (2itai-periodic) meromorphic continuation 3^™ of the function 

belong to the set V = A U {— Z + + 2-7raiZ} (f/iey are double for $l(z) £ — N = 
-1-Z+). 

b) If$lk > 1/2, ZTierc 

(1.25) lim (a/4) fe-1/,2 3* m (£;) = = r(fc + l/2)C(As + 1/2), 

(1.26) . lim a fc-1 / 2 3i m (fc) = v^r(fc + l)C(ife + 1/2). 
For aZZ A;, 

(1.27) (^) fe ~ 1/2 3g n ( fc ) = (4a) 1/2 - fe tan(7rfc)r(fc) 2 +3(A;), 

OS fl — ► CO. 

c) The function 3q m (Zc) tends to a nonzero limit as -R/c — > oo. Given 
k G — R+ stic/i i/tai {k + iR + Z} n "P = 0, i/te Zirrai 

(1.28) i/j(v;k) = lim 3g n (« — m + iv) 
existe and zs a continuous 2ira-periodic function of v £ R. 

The term 3 in (1-27) can be skipped for !RA; strictly less than 1/2 due to 
the a-factor. However, the formula with this term gives more than was stated 
in the theorem. In the critical strip { — 1/2 < ^tk < 1/2}, the right-hand side 
approximates well the left-hand side even when a is not too big . For instance, 
(a/4) fc ~ 1//2 3^ n is getting very close to 3 when *zsk ~ 10 or more to make T(/c) 2 
small enough. It is the same phenomenon (with the same T-limit) as in (1.20). 
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The set of admissible k in c) is dense everywhere in — R+ because so is 
V + iH + Z C C (for generic a). However the function ip{v] n) is well-defined. 
It is related to the existence of the totally singular zeta Q q considered above. 
In fact, ip is a limit of a certain sum similar to (1.3). 

Thanks to c), the number of zeros of 3q U (k) in the vertical strip 

{k — m — 1 < KA; < k — m} modulo 27raiZ 

equals the number of poles of 3^™ {k) (counted with multiplicities) in this strip 
for sufficiently big m. It is likely that even the positions of individual zeros 
(modulo Z + 2-7raiZ) approach certain limits when m — ► oo. This can be seen 
numerically. 

The number of zeros of 3^ n in the half-planes Kfc > C for any C G R 
is finite, which readily results from a). It is not difficult to calculate ip(iy;K) 
numerically for a ~ 10 or less. The change of aig(ip) / (2n) from v = through 
v = 2ira combined with the periodicity of 3g n can be used to calculate the 
exact number of its zeros modulo the period 2irai in the half-planes $lk > C. 

We can ensure the convergence to Riemann's zeta for all k switching to the 
"plus-zeta functions" defined for (exp(— x 2 /a) + l) -1 . In the classical theory, 
this substitution adds a simple factor to (, terminating the pole at s = k+1/2 = 
1 and improving the analytic properties. The zeros for -Rfc / 1/2 remain the 
same. Setting 

(1.29) 3!S(ife) d = H) / (exp(-x 2 /a) + iy 1 5 k (x;q)dx for m > 0, 

Jo 

(a/4) fe - 1 /2 3 ;-(A ; ) _ 3+ (fc) = f (l-2 1 / 2 - k )3(k)- 

The limit holds for all k G C except for the poles if 3!^ is replaced by the 
meromorphic continuation 3+^(&)- 

Locating the zeros. The numerical experiments demonstrate that 3g m 
has only two zeros ( mod 2irai) in the right half-plane $tk > e at least for 
e = 0.05, a < 25(q < 0.96). These zeros are connected to each other: 

ZQ = £ + Vh z\ = £ + (2™ - r])i, 7] < na. 

In a sense, these zeros appear because of the poles of 3g m at k = 0, irai. They 
vanish for 3!^(&). 

The zo slowly approaches /c = 1/2, the pole of C(k + 1/2): 

z (a = 3) = 0.0215 + 1.2746i, z (100) = 0.4419 + 0.7216i 

(1.30) zo(1000) = 0.4741 + 0.5560i, z (2000) = 0.4788 + 0.52081 

Switching to 3+^, we hope that its zeros in the strips > — e, e > 
are always deformations of the zeros of ((k + 1/2) for sufficiently big a > a(e). 
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This might be true for 3g™ as well with a reservation about the "low" zeros 
going towards k = 1/2 (similar to (1.30)) and their reflections. 

Let us mention without going into detail, that the following "theta"- 
deformation 

n QQl OO 

(1.31) 3j(fc) d = (-i) / CS2exp(-(nx) 2 /a))5 k (x;q)dx for $tk > 

J ° n =l 

also has only low zeros in the same range of a. Here the functional equation 
for the theta is used to integrate near 0. The analytic continuation of 3^ to 
the left half-plane is not known. 

We note that the roles of the left and right half-planes are not symmetric 
in the ^-theory. For instance, the zeros of 3+^ mainly appear in the left one. 
There is an explanation for this. 

Given a zero k = z £ iH of C,{k + 1/2), the linear approximations z+(a) 
for the zero z + (a) of 3+^ associated to a given zero k = z of Q(k + 1/2) can be 
calculated as follows: 

4(z + 1/2)C+Qz + 3/2) -{z- lX+jz - 1/2) 
Z+{a) ~ Z{1 + 12aC'(* + l/2)(l-2V2-*) } 

(1.32) for C'(«) = d((s)/ds, ( + (s) = (1 - 2 1 - s )((s). 

The formula results from (1.29). However we did not check the necessary 
estimates, so it is not completely justified at the moment. The first z+ (a) 
which we found in the right half-plane corresponds to z = 1977.2714i, well be- 
yond the range where we can calculate 3+^ numerically. After this point, the 
number of linear deformations in the right half-plane grows but not quickly. 
Generally speaking, they appear when the distances between the correspond- 
ing consecutive zeros of ( are getting too close, i.e. exactly when the linear 
approximation could not be reliable (numerically). 

The same strong tendency holds for the "real" zeta 3g £ - However now the 
deformations of the zeros prefer the left half-plane. The formula for the linear 
approximation z(a) of the zero z(a) of 3 r q e {k) corresponding to a given root z 
of C(k + 1/2) reads: 

U33) z(a)-z(l 4(z + l/2)C(. + 3/2)-(,-l)C(.-l/2) 
[L66) Z[a) ~ Z[l 12aC'(z + l/2) 1 

Note the opposite sign of the term after 1 and the missing factor (1 — 2 1 ~ s ). 
The formula for the plus-counterpart 3+ 9 (fc) of 3 r q e (k) is (1.32) with the minus 
after 1. The tendency to move to the right is practically the same for 3+ ? (&)- 
The factor (1 — 2 1_s ) does not contribute too much. 

Formula (1.33) gives a very good approximation for the zeros 2" (a) of 3 q (k) 
as well. Of course a must be big but not enormously big, because eventually 
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3>q(k) will go to its T-limit, which has no zeros in the critical strip. For instance, 
1(100) = 0.7091 + 14.1326* for z = 14.1347*, the first zero of ((k + 1/2). The 
exact zero zN(100) equals 0.7041 + 14.07791 Enlarging a: 

z(1000) = 0.0709 + 14.1345*, z tt (1000) = 0.0709 + 14.13401 

The first zero z of £(fc + l/2) with the "real" linear approximation zin the 
left half-plane (which we found) is z = 1267.57061 The next one is 1379.68331 
These zeros are beyond the range we are able to analyze numerically, although 
the calculations with 2>\ are 10 • • • 100 times faster than those for y q e , 

The a-evolution of the jj-zeros is very nontrivial. It is not too difficult to 
trace the evolution of the first ones back in a (at least, through a ~ 1). For 
instance, if z = 14.1347* 

(1.34) z fl (35) = 1.92 + 13.71*, z tt (ll) = 4.25 + 11.36*, 
z*(5) =0 + 8.72i, z\2) =-2.98 + 5.011 

Concerning big a, we were not able to reach the moment when the zero un- 
der consideration leaves the vicinity of z = 14.13471 The linear approximation 
works reasonably well at least untill a = 50000. One may expect to go to the 
vertical line Jftfc = 1/2. Employing the formula (1.20), ~ 0.0015 + 14.1320* 
when log(a) = 122.25, so it is still close to the classical zero for such a. When 
log (a) = 1020.34, approximately z$ ~ 0.4361 + 14.13341 This calculation 
is of course qualitative (as well as (1.20)), but presumably demonstrates the 
tendency. 

There are several (more than 1) examples of "low" zeros of 2>\ traveling 
from the left half-plane towards the natural "attractor" k = 1/2. The evolution 
of the first of them (a counterpart of the zero from (1.30)) is as follows: 

(1.35) a = 5 -J = 0.037 + 0.735*, a = 2733 : z tt = 0.409 + 0.477i, 
a = 4500 -J = 0.418 + 0.460i, a = 15187 : z^ = 0.436 + 0.424i, 
a = 76882 -J = 0.452 + 0.383i, a = 172984 : z» = 0.459 + 0.3651 

The convergence is very slow. It explains to a certain extent why the defor- 
mations of the zeros of £ stay put. The speed of decay at big a is expected to 
be much smaller for the "high" zeros. 

Summarizing, the distribution of the zeros of the imaginary g-zeta func- 
tions looks surprisingly regular. As to the 3g and its plus-counterpart, the 
main objects of our computer simulation, the g-zeros are also far from random 
in the horizontal strip {— 2e a < Qk < 2e a } for e a = ^pna (see (1.14)). 

There must be a reason for this. Numerically, a little decrease of the 
accuracy in the formula for 3g (say, ~ (5 • • • 10) ^/a terms in the range k = 
14i • • • 50i instead of the necessary ~ (15 • • • 20)y/a for a ~ 500 • • • 1000) and 
one can see how irregular the zeros could be. 
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2. Exact formulas 

First we will drop the factor exp((/c + 1) 2 /(4a)) — 1) in the formula (1.21) 
for 3q(k) (see (1.12)) and calculate the resulting sum. It is a variant of the 
celebrated constant term conjecture [AI,M1,M2,C3]. The latter corresponds 
to the imaginary integration over the period. The following theorem readily 
results from the (/-binomial theorem. Let e > 0, 

1 f OO — £2 

(2.1) € t q{k) d g 5 k (x;q)dx = 

a J oo+ei 



avr (1 - e -^'+ fc )/ a )(l - e -(^ fe )/ a ) 

~ T 11 (1 - e -(j+2fe)/a)n _ e ~(j+l)/a) X 
j=0 v ' v ' 

o^^ (1 _ e -a + 2 fc -l)/a )(1 _ e (/ +fc )/ a) 

{ ■ > Z^ll (l_ e -(l+k-l)/a)n_ e l/a) ■ 

j=0 1=1 v / v / 

The last series is convergent for all A; except for/cG— 1/2 — Z + + 27raiZ. The 
coincidence of the integral and the series holds provided that k/2 is inside the 
contour of integration. We will always use the series if k is outside. 

Theorem 2.1. For all k such that $lk < 0, - -^£q(k) = 

o fi (l- e -0-+fc)/ a )(l- e -0--fc)/ a) ^ a ^ l- e -(^)/a _ 
Z ii _ e -(j+M)/a)(l - e -0'+l)/«) 4^ ^ l-e-'/ a 

~ (1 _ e -(j'+fe)/<»)(l - e -(j+fe+l)/a) 

( 2-3 ) 2 11 (i-e-O'+afcV^n.e-a+iVo) • 

TTie /irm£ o/ i/te Zatter expression as a — > oo is p^r^fc+i) ■ 

Integrating the Gaussians. Now let us integrate the Gaussians g ±:c2 
with respect to 5 k - See [M1,0,J] about the classical (i.e. differential) multi- 
dimensional generalizations. The formula for the imaginary integral can be 
deduced from [CI] (in the case of A\). Let us discuss it first. 

Theorem 2.2. Provided that Sft/c > 0, 

i r+ooi 00 1 _ p -(j+k)/a 

(2.4) &™(k) d ^-J^ e x l a 5 k (x-q)dx = ^\[- i ' 



j=0 



e -(j+2k)/a ' 



(2.5) a k - ll2 &™{k) -► v^^^ as a oo. 
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Actually the theorem can be considered as a variant of the Ramanujan 
integrals, namely (3.12) and (3.13) from Askey's paper in [A]. They correspond 
to the two special cases k = ±oo, when 5 k has either no denominator or numer- 
ator, but involve more parameters. Richard Askey writes, "try to imagine what 
someone with Ramanujan's ability could do now with the computer algebra 
systems Of course the computers are useful for algebraic (/-calculations, 
but still supporting the (/-analytic theory (meromorphic continuations, zeros 
etc.) seems a more promising application. The theory is not mature yet and 
we need numerical experiments. However, first, computers must become at 
least 100 times faster. 

We note that the integral from (2.4) does not give the analytic continuation 
of & t J n (k) for Rk < 0. There is a simple way to extend the theorem to 3t.k > 
— 1/2 replacing 5 k by its asymmetric variant 

fr (1 - qi+ 2x )(l - g 3 + 1 ~ 2a! ) 

( 2 - 6 ) Vk{x;q) - 11 ^ _ qj+k+2 x^l _ q j+k+l-2xy 

The //-measure has many applications and appears naturally in the proof of the 
Macdonald constant term conjecture (formulated in terms of 5 only). Actually 
it is not surprising because the symmetry x ^ 1/2 — x of fi is quite common 
for many (/-functions (including the theta functions). The connection with the 
(5-measure is through the symmetrization: 

/+ooi i I —k/a />+ooi 
/(x 2 V fc dx = — — / f(x 2 )5 k dx for any /. 
-ooi ^ J —ooi 



The function 



1/4+ooi 



(2.8) = f - — - / e* 2 /y k (x; q) dx 

1 + e- k / a Ji/A-ooi 

coincides with <8>* m (A;) for > and gives its meromorphic continuation to 
K/c > —1/2 with the poles {irai + 2-iraiZ} coming from the factor before the 
integral. Respectively, (2.4) holds for such k if &™(k) is replaced by <5q(k). 
The passage 5 — > fi will be also used later for the zeta. 

Real integration. A more interesting exact formula is for 

, , i roo — ei 

(2.9) &l(k) = f - / e~ x /a 5 k (x;q) dx = 

^ % Joo+ei 



I oo+ei 

oo 



GTTyr (1 - e -0'+ fc )/ a )(l - e -(i- fc )/«) 
3 



oo 



(2.10) E^-^n 



L (1 _ e -0'+2*)/a)(l - e -(i+i)A») 

1 _ e -{l+2k-l)/a\M _ e {l+k)/a\ 



(1 _ e -(l+k-l)/a\(l _ e l/a) 
j=0 1=1 y J\ J 
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Here the convention is the same as for (t\(k): we use the integral if and only if 
it coincides with series, i.e. when k/2 is inside the contour of integration. For 
such k, (2.10) is nothing else but an application of the Cauchy theorem. 

Theorem 2.3. Let Q = exp(— 8a7r 2 ). Then for all k except for the poles 

k(k + ~\ \ SOL 1 - P -(i+ k )/ a 
*»(*) = ^)^ex P (M^ )n i_^__ for 

3=0 

<f>(k) = sin(^fc) J] 1 _ Qj/2 x 
3=1 

(2.11) (1 - e 2 ^V)(l - e - 27rifc Q J ')(l - e^ ik Q 2j ~ l ){l - e'^Q 2 ^) 

\/TTa 3fc 2 fc _fcsTT 1 — q J I 2 

= ~— q 2 (9 5 -9"')II(rr^)2 x 

oo 

(2.12) [](1 - ^'- fc )(l - ^ +fc )(l + ^72-1/4+^/2) (1 + ^72-1/4-^/2^ 
J'=l 

is a ^ oo ; 

(2.13) 0(ife) -► sin(vrA:), (a/4) fc_:l / 2 0}(ife) -► sin(vrfc)r(A; + 1/2). 



Concerning the limits, the convergence of 0(/c) to sin(7rA;) is very fast. 
They practically coincide even for a ~ 1 because Q is getting very close to 
zero. The Q 1 ! 2 is dual to q 1 ^ = exp(— l/(4a)) in the sense of the functional 
equation for the theta function: 

OO _| oo 

(2.14) Y, Q f/2 = (m/^k))(k = 0) = —= Y Q j2/4 - 



'j=— oo v 3=— oo 

Here we first check (2.11) and then use the functional equation to go to (2.12). 
The latter is straightforward. Combining (2.10) with (2.12) we can eliminate 
Q. Indeed, dividing by the common factors and distributing exp( fc ^"^ ) 
properly one arrives at the (/-identity: 

3=0 ^ 1=1 H 

yj (1 - g^ 2 )(l ~ ql +k )(l + gi/2-l/4+fc/2 )(1 + g3 72-l/4-fc/2) 

l\ (W) 

It is nothing else but the formula (1.11) from [CI] for the root system 
Ai and the polynomial spherical representation. The simplest particular cases 
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k = 0,1,2,1/2 readily result from the product formula for the classical ip = 
J2™=-oo c l n ( see (2-14)). It is possible that this identity for k € Z + /2 follows 
from known formulas. Then one can extend it analytically to arbitrary k. 
However here we establish (2.15) via (2.11). 
The key observation is that 

(2.16) (f>(k + 1) = -4>{k), (f>(k + 2p) = q 3kp+3p2 c/)(k) for p d = f 2-nai. 

The second formula results directly from the definition, the first one requires 
the technique of the shift operators which will be discussed later. These sym- 
metries and the normalization (2.14) are sufficient to reconstruct (p uniquely. 

A generalization. A natural setup for the theorem (in one dimension) 
is as follows. Let {p m (x)} for m = 0,1,2,... be the Rogers-Askey-Ismail 
polynomials (see [AI]) in terms of q nx for n € Z. They are x-even, pairwise 
orthogonal with respect to the pairing {/, g] = Const Term (fg5 k ) (where 
5k is replaced by the corresponding Laurent series), and normalized by the 
condition 

(2.17) = q mx + q- mx + lower powers, except for pf ] = 1. 

The multi-dimensional generalization is due to Macdonald (see e.g. [M2]). In 
this definition, k ^ —1, —2, —3, . . . , — m + 1 modulo 27raiZ. For instance, 

(2.18) p? =f + q~ x , # = ^ + + ^-^it^ - 



Setting 

. i i-co — ei 

&l(k;m,n) = f - / p£\x)pW [x)e~* /a S k (x;q) dx = 

Z J oo+ei 

ajL ™ (1 _ e -q+A:)/a )(1 _ e -b--fc)/a ) 
1 ' 2 11 (! _ e -(j+2fe)/a)( 1 _ e -(j+l)/a) X 

V p t« ( fc + j ) fMf fc + j zflga! A (1 - e-(^V)(l - e('+*>/«) 

I 2 /Pn 2 11 ( 1 _ e -(l + k-l)/a\f 1 _ e l/a) ' 

jr =0 Z=l v J\ / 

we deduce from [CI] the identity: 

(2.20) *!(*;m,n) = g - ^^^"^ ^(^^(^^(fe). 

The formula also holds for the (/-spherical functions, generalizations of 
Pm\pn^ to arbitrary complex m,n, and in the multi-dimensional setting (see 

2 2/ 

ibid, and [AI]). Moreover, the Gaussian J = q x = e~ x > a can be replaced 
by (/i/)*7 for the g-Hermite polynomials {hi(x)} (the Rogers polynomials at 
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k = oo) and the formal conjugation (q)* = q^ 1 . This is possible (in one dimen- 
sion only) because (/ty)*7 are eigenfunctions of the shift operator (see [BI]). 
We can make / complex too switching to the g-Hermite spherical functions. So 
a natural setting for Theorem 2.3 is with three complex parameters {k,l,m}. 
Hopefully one can add more replacing 5^ by an arbitrary truncation of the 
theta-function and involving BC\. Also (2.20) can be extended to nonsymmet- 
ric Macdonald polynomials. 

We note that both (2.4) and (2.15) can be deduced from known identities 
(Macdonald, Andrews). The general one involving Macdonald polynomials (or 
spherical functions) seems new. 



3. Analytic continuation 

In this section we construct meromorphic continuations of the "plus- 
counteparts" of 2>\ and 3g m and give the formulas for the a-limits. We mainly 
use the Cauchy theorem. There is another approach based on the shift operator 
technique, which will be applied in the next section. 

Let us start with the "sharp" function: 

V (k\d)- f - 6k ^ q ^ dx for 

(3.1) e = e+(a,d) = mm{J — ,— }, d > 0, m > -1/2. 

y ZiCL Zi 

The integration is with respect to the path {C e _ — C+} from (1.10). Recall that 
it goes from oo + ei to oo — ei through the origin. The function is well-defined 
and has no singularities in the horizontal strip = {—2e < 9/c < +2e} for 
?R:k > -1/2. 

Theorem 3.1. The meromorphic continuation of 3+ q (k\d) to all k is 
given by the formula 

air (1 - e -k'+ fc )/°) (1 - e ~0'- fe )/ a ) 

x 



--=nf 

3=0 V 



(1 _ e -0'+2fc)/o)(l _ e -(j+l)/o) 

~ 1 - e -U+k)/a * I _ e -(I+2fc-l)/a e xp(fcj/q) 

(3 " 2) 2- l- e -k/a 11 1 _ e -l/a 

]=0 1=1 ex PV 4a ) ^ 1 

The set of its poles in the strip K.\_ is {—1/2 — Z + }. It also has zeros at all 
integral k. For any k £ C except for the poles, 

a— >oo 4 

(3.3) 3+(k\d) d = f d-^ 2 - k S(k){l - 2 1 / 2 ~ k ). 
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Formula (3.2) results from Cauchy's theorem. Thanks to the exp(dx 2 /a), 
the convergence is for all k (except for the poles). The zeros of 3 +q (k) at 
k £ Z_|_ come from the product before the summation. It is less obvious for 
negative integers (use Theorem 4.2 below). Actually it is a general property 
of the difference Fourier transform (cf. (2.20)). If $tk > —1/2, the limit (3.3) 
is due to the integral representation. For the other k, one can use the shift 
operators (the next section). We note that the calculation of the limit for 
$lk < —1/2 is equivalent to that for the imaginary q-zeta. 

Imaginary zeta. We define it for -Rfc > as follows (cf. (1.29)): 

(3.4) &™{k\d) = 3 +q {k\d) = f (-») f lQ °(eM-dx 2 /a) + ir 1 5 k (x;q)dx. 

Jo 



We will drop {im} till the end of the section. Its limit as a — > oo can be readily 
calculated: 

(3.5) (a/4) fe -V2 3+(?(A; ) ^ 3+(fc|d) _ 



The first step is to extend this function to small negative $lk. Let 
(3.6) e = min-L/^, -}, -2e <$lk < 2e. 



The following function is analytic for such k except for the poles at {irai + 
2iraiZ} and coincides with 3 +q (k\d) when $tk > 0: 

d+ q (k\d) =' - / (exp(-dx 2 /a)+l)- 1 5 k (x;q)dx + f (k), 

def an ^ (1 - e"^ +fc )/ a )(l - e"(^ fe )/ a ) 
(3-7) il(k\d) = —[[ - _ e _ a+2fc)/a)(1 _ e _ - + i)/a) x 

" (l- e -q+2fe-l)/a)(l_ e (f+fc)/a) ±2^ 1 

11 (I _ e -(l+k-l)/a)n _ e l/a) 2^ / -djk+n+2-naim) 2 n . -■ " 

/=! v l\ J m=-oo c-^-PV 4 a J 1 



To see this, we replace (— i) J* °° l in (3.4) by (— i/2) and then by 

(—i/2) f**^+$o provided that < $lk < 2e. Here we applied Cauchy's 
theorem 
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Theorem 3.2. Let £, = ((21 + l)vrai/d) 1/2 , where 1 = 0,1,,... 

S d = n U A U A+, n = {-Z+/2 + naiZ} \ {-Z+ + 27razZ}, 
A = {-2^ - Z+ + 2iroiZ} U {-2^f - Z+ + 27raiZ}, 
(3.8) A + = {2& - Z+ + 27raiZ} U {2ff - Z+ + 2ttcmZ}. 

5y * we mean the complex conjugation. The series 

oo oo 

y +q (k\d) = J2*n(k\d), 3l q (k\d) = Y,2l(k\d), 

n=0 1=0 

S(k\d) = a7T (r 1 TT (i-e-^' )/B )(i-e-^)/") 

* l { 1 ' ' " 2d V* 1 Al (i _ e -(j+fc+2^)/a)(i _ g-y+fc-^J/a) 



^ ^ li ^ _ e -0"+fc+2£f )/«)(! _ e -(i+fe-2€D/a) j 

are absolutely convergent for any fixed k S with Jftfc < 0. TTie difference 

(3.10) 3^ fl (fc|d) "= 3%(fc|d) - d x +q m 

represents the meromorphic continuation of 3 +q (k\d) to all 9ftfc < except for 
the set of the poles V = II U A. 

It is not difficult to check that the poles of 3+q and 2>\ q apart from V have 
the same residues and cancel each other in (3.10). The latter set is discrete 
in contrast to S which is dense everywhere (for generic a). We will skip the 
convergence estimates because the theorem can be reproved using the shift 
operators without termwise consideration. 

Actually the self-termination of the poles readily follows from the integral 
representation (3.7). The latter coincides with and has no poles for > 
Sftfc > — 2e. If the differences of residues of the poles from A + are nonzero they 
can be seen in this strip (for generic a). This argument and the theorem can 
be somewhat generalized as follows. 

Proposition 3.3. Given r € Z + , the meromorphic continuation of the 
integral 3+^ (see (3.4)) from the strip K r = {k, —r<^Stk<—r + 1} to all 
negative $ik is given by the series 

oo 

(3.11) 3? q {k\d) = f Yl s ^n(j-r + 0.1) ^(k\d) - 3 x +q (k\d). 

3=0 



Its poles QRk < 0, r > 0) belong to V r = {-r + A} U II \ {1/2 - r + iraiZ}. 
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This statement will be involved to analyze 3+g 3 as — > — oo upon the 
following renormalization: 

(3.12) 3+ q (k\d) = P{k)3 +q {k\d), 3+,{k\d) = p(k)3Z,{h\d) for 

( 3.13) «^n, (1 -C"'^!v 

These functions also have meromorphic continuations to £ C. The additional 
factor simplifies the singularties. The poles of the continuation of 3+<j (all k) 
form the set V = A U { — 1 — Z + + 27raiZ}. The poles of the second function 
(negative Wt, r > 0) constitute V r = {-r + A} U {-1 - Z + + 2iraiZ}. It also 
has "two" zeros {1/2 — r + iraiZ} modulo 2iraiZ. 

Similar methods can be applied to the 3q m from Section 1 corresponding to 
the standard kernel (exp(— x 2 /a) — l) -1 . One has to replace $„(k\d) switching 
to the factors (exp(— (k + n + 27raim) 2 /(4a)) — 1) and use & = (lirai/ 'd) l l 2 
starting the summation from 1 = 1. Other ingredients remain the same. 

k-limits. The behavior of g-zeta functions as k — > oo is very non-classical. 
Let us first tend to +oo. We may use the integral representation: 



(3.14) 3 +q (k\d) - V+ = H) / 

■/ o 



nr =0 (i-^' +2;c )(i-^- 2;c )^ 



exp(— dx 2 /a) + 1 

as ?R.k — > +oo. So the limit is a positive constant. 

The calculation is more involved when $tk goes to — oo. Namely, the 
integral representation cannot be used anymore and the renormalization (3.12) 
is necessary. We will apply Proposition 3.3 as r — > oo and k = k — r + vi for 
< n < 1. First of all, 

(3.15) «) <- 3;j(fc|d) - -3+,(fc|d) - 3i,(fc|d) 

for a continuous function </>(zv; k) which is 27ra-periodic with respect to v € R. 
Here we use the integral representation for 3+q(k\d) and evaluate the change 
of the integrand as k i— > A; — 1 (the multiplicator) for big — 3t.k. Thanks to the 
renormalization, the integral tends to a Z-periodic function of k. Of course it 
is necessary to avoid the singularities k £ S. 

Thus we established a relation between 3+ g and 3+ g in the considered 
limit. Let us use it: 

(3.16) lim3^ = lim(3+ ? - 3+ g ) = -<f>{v; n) - 21im3+ ? (A;|d). 
If k is fixed (and is in a general position), 

(3-17) 2>Uk\d) = -[ ^ 5 f 2 \ q \ dx , 

v ; 2i y« exp(-dx 2 /a) + 1' 
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where the integration contour <C is the boundary of a neighborhood of the 
union 

{(6 + R + )(l + *)U(e + R+)(l-*)}. 

The orientation is standard. The neighborhood must contain no points from 
±{(k/2 + Z+/2) + iraiZ}. So it depends on k, v and approaches the diagonals 
for big $lx. Now we can go from k to k — 1 evaluating [3{k — l)<5fc_i(/3(&;)<5fc) -1 
as r — > oo. 

Theorem 3.4. Given < k < 1 suc/t t/iai {k + iR + Z} n P = 0, the 

limit 

(3.18) iP{v,k) = hm Z+a {k- — r + iv) 

exists and is a continuous 2ira-periodic function of v G R. 

The theorem is the "plus- variant" of Theorem 1.4, where V was somewhat 
different. The set of admissible k is dense everywhere in [0, 1] because so is 
V + iH + Z C C (for generic a). However the function k) is well-defined. 
It is directly related to the existence of the totally singular zeta Q from Section 
1, because 3+ g is nothing else but a variant of the function from (1.7). 

Thanks to the theorem, the asymptotic number of zeros N r of Z r ^q{k\d) 
for big (generic) r in the box 

def 

B r = {k — r — 1 < Jft/c < k — r} modulo 2iraiZ 

coincides with that for Z+q{k\d) and equals the number of the poles, i.e. the 
points from VC\B r . Indeed, it suffices to calculate the difference of the changes 
of axg(ip{v)) over the "periods" at k = and k = 1. Note that N r is quadratic 
in r : 



(3.19) N r = 2 



1 + d(r + 1 



\2 



4a7T 



+ 1, < k < 1, r > 1. 



It is likely that the positions of individual zeros (modulo Z + 2iraiZ) 
approach certain limits when r — > oo. To ensure this, the convergence to 
ip(y; n) must be fast enough. This can be seen numerically. For instance, if 
a = 1 = d the function 3 +q (k\d) has the following (unique?) stable real zero: 

z e = -6.668819069126, z 7 = -7.669834033011, 
(3.20) stable zero z r = -r - 0.67042591089415. 

The stabilization is rapid (for small a). At least one real zero must exist in B r 
for sufficiently big r because N r is odd. Recall that 3 +q (k\d) is real on R, so 
the complex zeros appear in pairs. 
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a-limits. Let us establish the relations to the classical counterparts of 
the considered zeta functions: 

3+(fc|d) = d- 1 / 2 - fc (l-2 1 / 2 - fc )3(A:), 

(3.21) 3+(fc|<Q = ^ 1/2 - fc (l-2 1 / 2 - fc )^r(A:+l)C(fe + l/2). 

See (3.3) and (1.26). 

Theorem 3.5. Let 3+q, 3+^ &e i/ie meromorphic continuations of 3+q 
and 3+q to a// fc constructed above. Then 

(3.22) lim (a/4) fc -V 2 3 + (k\d) = 3 + (k\d) for k ^ -1/2,-3/2, ... , 

(3.23) lim a fc ~ 1/2 3+ g (A:|d) = 3+(fc|d) for k + -1, -2, . . . . 



The proof is based on the integral representation (3.4) for > 0. As to 
3+ g , we also use (1.5) from Lemma 1.1. When 5ft/c > —1/2, we can go from 5k 
to Hk (see Section 2 and below) or involve 3+ g from (3.7). Let us deduce the 
formula for the limit of 3+^ for < — 1/2 from (3.3), Theorem 3.1. Generally 
speaking, the real and imaginary (/-zeta functions are not connected, because 
the contours of integration are different and the signs of the Gaussians are 
opposite. However, there is a relation in the limit. 

First, let us consider 3>\ q {k\d) (see (3.10) and (3.25)). Since a — > oo, the 
terms in the sums (3.7) for ^(k\d) with m > are very small and can be 
ignored. Therefore 



+ ^ V aVT-fi (1 - e -Q+*V») (1 - c - » , 

+ - 2 JUL _ e -(j+2fe)/a)(l _ e -0'+l)/a) 

n 



-(j+fe)/awi _ -(j-fc)/a\ 



i=0 i=0 

(1 _ e-C+afe-^/a^n _ e (J+fc)/a\ 1 



, , (l-e-('+*-i)/«)(l-eV«) exp( ^±2l! ) + i 

(3.24) -C g (fc)+3%. 

Here we used the formula for the "sharp constant term" from (2.2) and of 
course the identity (exp(x) + 1) _1 + (exp(— x) + 1) _1 = 1. 

Since <t\(k) = 0(a) (Theorem 2.1) and 3+ 9 = 0(a 1/2 ~ fc ), the latter will 
dominate for big a because $lk < —1/2. Thus 3+q 9 — 3+ g — 3+ 9 , and thanks 
to (3.3) it suffices to calculate the limit of 3+ g - This can be done via the Moak 
theorem (Lemma 1.1) exactly as we did in (1.6) with ( q : 

oo 

(3.25) Z x +q (k\d) - £ ^(a/4)- k ((ax l )- 1 ' 2 (-x l ) k + (-ox,)" 1 ^^)*) 
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for x\ = ({21 + l)ni/d). Then we use the classical functional equation. 

Summarizing, our technique of analytic continuation of the imaginary q- 
zeta is a q- variant of the so-called Riemann's first proof of the functional equa- 
tion (see [E], 1.6 and the appendix, the translation of Riemann's "On the 
number of primes less than a given magnitude" ) . In contrast to the classical 
theory, the resulting functions are periodic in the imaginary direction and have 
limits as — > ±oo. The "sharp" g-zeta functions 3? and 3+ g can be repre- 
sented as convergent series for all k, which has no counterpart for Riemann's 
zeta. 



4. Shift operator 

The shift operator in the setup of this paper is the so-called Askey- Wilson 
operator, the following ^-deformation of the differentiation: 

Sf(x)= f (q x - q- x )-\r - t" 1 )/^), 

(4.1) r/(x) = /(x-i), g = exp(-l/a). 
For instance, 

(4.2) Sj = -q lf S, Sj- 1 = g- 1/4 7 _1 for 7 = q x \ 

The name is "shift operator" because its action on the hypergeometric func- 
tions, both classical and basic (difference), results in the shift of the parameters 
[AW]. For example, 

(4.3) Spg\x;k) = r m/2 -« m/2 )p£i(^Hl), m = l,2,-.. 

for the Rogers polynomials from (2.17). 

The differential shift operators for arbitrary root systems are due to Op- 
dam (see [O] and also [H] for the interpretation via the Dunkl operators). The 
general difference ones are considered in [C3]. They depend on k in contrast to 
the simplest case considered here (Ai). The celebrated constant term conjec- 
ture and the Macdonald-Mehta conjecture were verified using these operators 
(see [Ml, O, C3, CI]). They were also used in [O] for analytic continuations. 

In this section, the shift operators will be applied to go from positive K/c 
to the left without involving Cauchy's theorem and convergence estimates of 
the pole decomposition. Some properties of the difference shift operators have 
no counterparts in the classical theory. For instance, they give a connection 
between the values of (/-zeta functions at k and k + 1 which collapses in the 
limit. We do not discuss the technique of the shift opeartors in detail. It is 
similar to the algebraic theory [C1,C3]. See also [C2]. 
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Shift-formula. Let na > e > and f be one of the following integrations: 

p j i-Kai/2 p re+coi 

I = : / dx, = T . dx, 

Jet ' Ka ' 1 J -nai/2 Jim ^ l J e—ooi 

p -i p-\-oo — ei p i roo — ei 

(4.4) = o dx, =- dx. 

J re ^ J—oo—ei J(j 2i J oo+ei 

Given it > 0, a function E^x) is called u-regular, if it is even, analytic in 
D ct = {x\ - u < "fox < u} for CT, 

(4.5) Di m = {x\ — u — e < 'Rx < u + e} for / , 

J im 

D re = {x G ±ei + R} for / (no u), 

J re 

Dft = D re U {x\ — e < < e, — u < $tx < u} for / , 

h 

and continuous on the closure of D. It is also assumed to be 7rai-periodic for 
CT and have coinciding continuous limits 

e ± (K + i'i) = lim E(x)5 K+vi (x) 

in the case of the imaginary integration f.. 

The following proposition results from (4.3). The proof is similar to [C2] 
(in the one-dimensional setup). 

PROPOSITION 4.1. Let K e be the set of k such that %lk > 2e for the 
imaginary integration, > for CT, {^kk > 0, \Qk\ < 2e} for the sharp one 
L, and K e = {^k + 2iraZ / ±2e} for the real integration. Given a 1/2-regular 

(k) 

E and a Rogers polynomial p m (m > 1), 

(4.6) J S{E)pl + _^5 k+1 = q k (q' k -f- q ^T) J E pV>8 h 

for k 6 K e , provided the existence of the integrals. 

Note that we do not take iR as the path for the imaginary integration 
because S(E) may have singularities at 2-7raiZ. Concerning the convergence 
of the integrals from (4.6), the function E can be arbitrary (say, continuous) 
for CT. For the imaginary integration, we can use the Stirling formula from 
Lemma 1.1. For instance, the integrability of |x| _25Rfc in the imaginary 

directions inside Di m is sufficient. The estimates in the real direction are more 
involved but not difficult either. 

Next, we put 



(4.7) 5 k+1 = (l-g 2 * +fc )(l-g- 2 *+ fe )5, 



k, 
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express (1 — q 2x+k ){\ — q~ 2x + k ) i n terms of (see (2.18)), and then apply 
the proposition twice to turn into a constant. 

Theorem 4.2. Given e > 0, let E(x) be 1-regular. Assuming (again) 
that k G K e and the integrals are well-defined, 

(l-g 2fc+1 )(l-g 2fc ) f pr 
r a 3/2+k r 

( 4 - 8 ) y es ™ + (1 _ g2fc ; 2)(1 _ g2fc+ 3) y ^)w 



The generalization of the theorem to arbitrary root systems is straight- 
forward. However the explicit formulas for 5k+i/5k hi terms of "small" Mac- 
donald's polynomials are getting complicated. 

An immediate application of the theorem is the constant term conjecture 
in the Ai-case (see [AI,M1] and the papers of Milne and Stanton from [A]). 
Let / = j cV E = 1. Then 



d k = t, — — 3ittt / d k 



ct 
oo 



(l-q 2k )(l-q 2k+l )Jct +1 



(49) 2n q-^)(i-3 fc+j+1 ) 

The resulting product formula holds for 3ft/c > 0. The reasoning is the same 
for other integrations with a reservation about the convergence conditions. For 
instance, we get (2.3) for the sharp integration. However now the convergence 
is for Ikk < 0, so we need to employ a variant of the theorem for the negative 
countours of integration and use the chain 

Actually this chain can be applied to establish (4.9) as well. First, we as- 
sume that k £ Z + , second, get the exact formula and, third, use the analytic 
continuation. 

It is not surprising that the products in (4.9) and (2.3) coincide. Formally, 

-nai \ Ebk = 2i E5k 
Jet J$ 

if E is regular everywhere. 

Another application is to the Gaussians. Using (4.2) we come to (2.4) for 
E = 7 _1 . When the integration is sharp and E = 7, we establish the symmetry 
(2.16), which played the key role in Theorem 2.3. 
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Let us apply S 2 to our real and imaginary plus-kernels E^ = (l+q ±x2 ) 1 
when d = 1 : 

(4 11) S 2 (E^) = ± - ^ 

So it is the same for the imaginary and real integration up to a sign. The 
complete factorization is a special feature of S(E^) and S 2 (E^). The for- 
mulas are not that nice for d ^ 1 and for higher S'-powers of E^ . We mention 
that big 5-powers have remarkable stabilization properties which we hope to 
discuss somewhere else. 

Analytic continuation. The theorem can be applied infinitely many 
times if the function E is oo-regular, for instance, in the case of the real inte- 
gration J re . Of course it is necessary to ensure the convergence of the involved 
integrals and the resulting series. It may lead to interesting transformations 
and identities. However it seems that the most promising applications are an- 
alytic. Namely, using (4.8), we can extend the function j E5 k to negative 
when E is analytic in the entire complex plane or proper vertical strips. If E is 
meromorphic the method can still be used. Of course the poles will contribute 
to (4.8), but the extra terms are meromorphic for reasonable E. 

Corollary 4.3. In the setup of Theorem 4.2, let (£(k) d = J E(x)5 k for 
k £ K e . The function 



(l- g 2k+l)(l_ g 2fc) 

r S 2 (E)S k+2 



q 3/2 +k{1 _ qk+ 2 ){1 _ qk+ l ) 



(1 - g 2fc + 4 )(l - q 2k + 3 )(l - q 2k + 2 )(l - q 2 ^ 1 ) 



is a meromorphic continuation of <£(k) to K e — 1, with f S 2 {E)5k+2 being 
analytic. 

Assuming that E is r-regular for r G 1 + Z + and all the integrals are 
convergent, the iterations of (4-12) provide a meromorphic continuation of 
<E(k) to K e — r with simple poles belonging to the set 

{-Z+/2 + iraiZ} \ {-Z + + 2vraiZ}. 



The described procedure becomes somewhat sharper for the imaginary 
integration when 5 k is replaced by ji^ from (2.6). Namely, 



(4.13) 



C(fc) = 2(l+g fe )- 1 / E(x)n k 

J im 
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for k £ K € , where the latter integral gives a meromorphic continuation of <£(k) 
to K e — 1/2 with the poles in nai + 27raiZ. Respectively, substituting /j, for S, 
the corollary provides the continuation to K e — r — 1/2 for r-regular E. 

Note that the above set of poles coincides with n from (3.8). Thus the 
results of the previous section can be reproved by means of the shift operator 
technique for E(x) = (exp(— dx 2 /a) + 1) _1 with a reservation about meromor- 
phic E. Using (4.12) only, we have to stop once the first pole of E appears, i.e. 
when m = -25ft(£ ) (Theorem 3.2). 

Limits. An important application of the shift operator is to the fc-limits 
and a-limits. We will discuss here only the latter. Formula (3.3) is still not 
justified for k < —1/2. Now we can proceed as follows. Let us take E = 
(exp(dx 2 /a) + l)- 1 . 

The diagonal poles are departing from the origin as a — ► oo, so we can 
define the analytic continuation using <£, <B^\ <£ ( - 2 - ) etc. (any number of them). 
It suffices to consider the first step, the passage from 3ft/c > —1/2 to 3ft/c > —3/2. 
Then we can repeat the procedure. For big a, 

(£)*-iA e c)( fc ) =* Jt±lL. { ^ + m e(k + l) + 

< 4 ' 14 > (* + 3/2)Wl/») ( f )W 'i S ' (g) ^- 

We already know (thanks to the initial integral representation) that 

(lf+i/2^ k + !) _ s in(7rA;)3+(A; + l\d) for k > -3/2. 

Actually the exact formula does not matter since this term comes with the 
factor (1/a) and will vanish. So we need to evaluate the second term only. For 
the sake of simplicity, let d = 1. Using (4.11), 

lim^^/V 1 ^) = 

a^oo 4 



1 



2i(fc + 3/2)(fc + l/2) 



pOO — El 

/ z-V 2 (- z ) k+2 £Wdz, 

J oo+ei 



(4.15) = - e " (1 



(1 + e 2 ) 3 ' 



Here we substituted z = x 2 /a. The integration path must be of course 
recalculated. However it can be deformed to the same "sharp" shape for a 
small e > because k + 2 > 1/2. In fact, k + 2 > —1/2 is sufficient. Then we 
may integrate by parts: 8^ is the second derivative of S = (1 + e z )^ 1 . The 
latter is not surprising because S 2 is a deformation of the second derivative. 
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Finally, 

lim^^/V 1 ^) = - \~y\-zfSdz = 

(4.16) (1 - 2 1/2 ~ k ) sm(nk)C(k + l/2)r(ife + 1/2) = sin(7rifc)3 + (ife|l). 

The last equality holds for $tk > —1/2. The integral from (4.15) represents 
3+(fe|l) for $lk > —5/2 (we need it only when $lk > —3/2). So, indeed, we 
moved by —1 to the left. Similarly, we can use and so on. This proves 
(3.3). One can (try to) apply this method for finding other terms of the (1/a)- 
expansion of ^ +q {k\d) and other (/-zeta functions. Moreover, it may work as 
a — > and q tends to roots of unity, however this is beyond the framework of 
the paper. 

Summarizing, the technique of analytic continuation from this section has 
the following classical origin. Each integration by parts shifts by —1 the k- 
domain of the integral representation (4.16) for the "plus-zeta". There are 
of course other similar examples in the classical theory of zeta. However the 
difference integrating by parts, namely (4.12), seems more promising. It has 
better analytic properties and, what is important, establishes a certain con- 
nection between the values of (/-zeta functions (or more general integrals) at k 
and k + 1, collapsing in the limit. The stablity of our zeta functions modulo 
Z, which has no classical counterpart, is also directly related to (4.12). 



5. Locating the zeros 

In this section we discuss (mainly numerically) the zeros of (/-zeta functions 
and demonstrate that their behavior is quite regular at least for SR/c > —1/2. 
Almost certainly the evolution of the (/-zeros is related to the distribution of 
the classical ones. We will give some evidence. 

Imaginary zeta. Let us begin with the "plus-imaginary" 3+^(k\d) from 
(3.4). Recall that 

pooi -i 

(5.1) 3T q (k\d) = H) J e _ dx2/a + i S k (x;q)dx for &k > 0, 

lim {a/A) k - l l 2 #™(k\d) = 3 + (k\d) = (1 - 2 1 / 2 - fe )(T 1 / 2 - fc 3(A:), 

where 3(fc) = ((k + l/2)r(fc + 1/2). 

Its analytic continuation has the same limit for all k except for the poles. Since 
3!S is 27roi-periodic and has the limit (a positive constant) as $lk — > oo the 
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number of its zeros N+(a; d) in the strip = {k | 3ft/c > e > 0} modulo 2irai 
is as follows: 

(5-2) N%(a-d) = iarg(3V m g )r fc tr- 



Here we used that 3+^ is real on R + and nai + R+. 

Numerically, it has no zeros in the right half-plane at least for a = 25, e = 
0.05, d = 2. Here is the table of the arguments a(n) = arg(3!^)(7rain/21 + e) 
for n = 1, . . . , 20 (a(0) = = a(21)): 

(5.3) a(l - 5) = -7.45845 - 11.6586 - 12.1104 - 12.4365 - 11.5873 
a( 6 - 10) = -10.1478 - 8.26116 - 6.03690 - 3.56672 - 1.01060 
a(ll - 15) = 1.61364 4.19644 6.66558 8.91015 10.8286 
a(16 - 20) = 12.2974 13.1573 13.1854 12.0256 8.96469. 

The total number of points on the half-period is ~ 320 in this calculation. 
The behavior of the argument becomes more complicated approximately after 
a = 30. Our computer program (it adds points automatically to ensure the 
required accuracy) was not able to reach a reliable answer for such a. Eventually 
the zeros will appear in the right half-plane. Indeed, the g-deformations of the 
classical zeros with $lk = 1/2 coming from the factor (1 — 2 1 / 2 ~ fc ) must arrive. 

Without going into detail, let us mention that similar results hold when 
(exp+1)" 1 is replaced by regular modular functions (cf. (1.31)). The "imag- 
inary" g-deformations of the corresponding L-functions (with proper gamma- 
factors) are well defined for ^Stk > 0. Sometimes elementary functions can be 
used instead of modular ones in the integral representations. Say, for the L- 
function associated with the "standard" Dirichlet character modulo 3, we can 
proceed as follows: 

rooi -i 

(5.4) 3S(*I<*) = H) I e - dx y a + l + edx y a ^l)dx for Xk > 0. 



Its analytic continuation approaches the classical L-function up to a gamma- 
factor for all k 6 C as a — ► oo. It has no g-zeros in the right half-plane at least 
when a ~ 25, d = 2. 

Recall that in the critical strip 1/2 > ?R.k > —1/2, we use the analytic 
continuation, namely, either 

3 b +q = 3+/ from (3.7), or 

/•1/4+ooi 

(5.5) y +q = -i(l + q k )- 1 (exp(-dx 2 /a) + l)- 1 fi k (x;q)dx. 

Jl/4—ooi 
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They coincide (see (3.10) and (4.13)). The passage to negative is necessary 
because the most interesting zeros, the ^-deformations of the classical ones, 
prefer the left half-plane. The theoretical explanation is as follows. 

Given a zero k = z of the classical £(/c + 1/2) (only z £ zR+ will be con- 
sidered) we can calculate the linear approximation 2+ (a) for the corresponding 
exact zero z + (a) of 3+ g : 



We expand the integral representation in terms of (l/o) using a refined variant 
of Lemma 1.1. Another approach is based on the shift operator technique. 
However we did not check all the estimates in (5.6) so the exact range (both 
a, k) where it can be used is not determined at the moment. 

The first 2+ (a) which we found in the right half-plane corresponds to 
z = 1977.27141 The formula predicts that this z could have the g-deformation 
z + (a) to the right, but we did not even try to find it numerically. It is well 
beyond the capacity of our computer program and the speed of existing (avail- 
able) computers. The applicability of the linear formula for so big z is also 
unclear. A more systematic analysis of the Moak-Stirling formula is necessary. 

The negativity of 9ft (z+ (a)) is a certain property of Riemann's zeta and 
its derivative. As we will see, it is connected with known facts. However there 
is no explanation why the tendency is that strong. 

Sharp zeta. Let us switch entirely to the main object of our computer 
simulation, the plus-sharp (/-zeta function ^ +q {k\d) from (3.2). It tends to 
(a/4) 1//2 ~ fc sin(7r£:)3+(fc|cf) for all k except for the poles (Theorem 3.1). We 
will analyze it in the strip 



where it has poles at {—1/2 — Z + } and vanishes at all integral k. We assume 
that a > 2/(i:d 2 ) which provides that the poles with Q/c = ira will not appear 
in this strip. 

The formula for the linear approximation z+(a) for the exact zero z^_(a) 
of 3+o corresponding to a given zero z of ((k + 1/2) reads 



ones z in the strip K\ and the values of the corresponding linear approxima- 




) 



(5.7) 
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tions z\. Let d = 2, a = 750 : 



2=14.1347*' 4 


= 0.1304 + 14.1450* 


z\ = 0.1302 + 14.1465* 


21.0220* 


0.3514 + 21.0702* 


0.3504 + 21.0771* 


25.0109* 


0.5641 + 24.9586* 


0.5745 + 24.9643* 


30.4249* 


0.9046 + 30.4014* 


0.9134 + 30.4077i 


32.9351*' 


1.1051 + 33.0341* 


1.0998 + 33.0854* 


37.5862* 


1.6449 + 37.9660* 


1.7675 + 38.1895* 


40.9187*' 


1.9080 + 40.8119*' 


1.9141 + 40.7816* 


43.3271*' 


2.2860 + 43.2485* 


2.4497 + 43.3138i 


48.0052*' 


2.9259 + 47.8424*' 


3.1103 + 47.5578* 



We did not try to reach high accuracy. However the existence of the zeros 
was carefully justified by means of sequences of integrations over diminishing 
rectangles around the consecutive approximations. The same integrations were 
used to calculate the g-zeros themselves. 

Note that the standard Newton method and its modifications work for 
the first 5 zeros but diverge for the other 4 even when applied to the zeros we 
found via the contour integrating (z+). 

For a till 10000, there is no problem with finding the first zeros, but 
calculating the last ones is difficult. 

The zeros exist for small a as well. If d = 2 the least a with the zero in 
the strip is a = 61 (anyway a = 60 is smaller than necessary). The zero 
equals z\ = 1.1721 + 13.8088* whereas 2e = 13.8433. The corresponding linear 
approximation = 1.60159 + 14.2789* is beyond the if -strip. 

This does not mean that we cannot consider the (j-deformations of the 
classical zeros for small a. For instance, tracing 1.1721 + 13.8088* back in a : 

a = 30 : z\ = -3.6951 + 10.5620*', z\ = 3.2566 + 14.4279*', 
(5.10) a = 1.98952 : z\ = -11.5000 + 2.5004*'. 

The first is very far from its linear approximation z\ , the second seems to have 
nothing to do with the classical z = 14.1347* at all. However both are indeed 
the (/-deformations of the latter. 

We note that outside K+, the appearance of diagonal poles (from A) 
changes the picture of zeros (and their number) dramatically. For instance, 
2e = 9.70813 for a = 30, d = 2 and there is a zero 0.3862 + 9.9175*'. It is also 
necessary to take into consideration the classical zeros due to the factor (1 — 
2 i/2-fc)_ For instance) 0.8223+9.0518*', 1.9601+17.6542* are exact deformations 
of 0.5 + 9.0647i, 0.5 + 18.1294*' for d = 2, a = 100. 
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According to our calculations, it is possible that all zeros z + inside 
(or a bit smaller strip) and with SR/c > — 1/2 (a > 2 /(nd 2 )) are g-deformations 
of the classical ones ($tk = 0, 1/2). However this must be checked more sys- 
tematically. Vice versa, all classical zeros presumably have the deformations 
for all q, but maybe far from the critical strip (see (5.9), (5.10)). 

Why do the zeros move to the right? The last topic we are going to 
discuss is a qualitative analysis of the positivity of §ft(S"). A similar question 
is about ups and downs of the linear deformations. They look random in 
the table (5.9), but there are interesting tendencies, especially for the usual 
"minus-zeta" . We will not touch it upon here. 

Since z G iR+, we need to examine the positivity of 



(5.11) V+ (z\d) = f 



Ad- l {z + ±)(1 - 2- 1 /2-*)£( z + |) _ d ( z _ x)(! _ 2 3/2-*) C ( z _ 1 



)• 



V 12aC , (^ + l/2)(l-2 1 /2-^) 

If d = 2 the smallest zero z when (5.11) becomes negative is z = 1977. 27i. 
The next one is 2254. 56i. Actually this property has almost nothing to do with 
the zeros of £. Indeed, the first interval where i] + {z\2) takes negative values 
(Ssz > 6) is [1977.24i, 1977.35i]. Since there are many zeros of Riemann's zeta, 
there is nothing unusual that one of them appeared in this interval. The 
tendency is practcally the same from d = 0.2 till d = oo. Only around d = 0.1, 
a (much) smaller "negative" zero appears: 163.031 The next one is 353.49i. 

Thus the positivity phenomenon is mainly because of the coefficient of d 
in (5.11): 

(5 12) / -(z-l)(l-23/^)C( z -l/2) x 
[b ' ^\ 12a<'(z + l/2)(l-2V2-*) J' 

Here the factor (1 — 2 3 / 2 ~ z )/(l — 2 1 / 2-2 ) is not too significant. Numerically, 
it is because its real part, always positive, provides the major contribution. Let 
us drop it (together with 12a). So we need to interpret the tendency 

(5.13) V (z) = f - z)((z - l/2)/C'(* + 1/2)) > 

on the zeros z G iR+ of (. 

Note that r] + (z\d) without the factors (1 — 2"') describes the linear defor- 
mations in the case of the usual zeta (see (1.33)). Therefore a little influence 
of the skipped factors could be expected a priori. Actually we can simlify rj 
even more replacing (z — 1) by i, since it almost belongs to H + i. However this 
factor is necessary to fix the positivity for small zeros, so we will not touch it. 

The phenomenon is again not about the zeros. The first interval where 
rj(z) takes negative values is [1267.47z, 1267.70i] whereas the first zero with 



ON Q-ANALOGUES OF RIEMANN'S ZETA 



33 



negative rj is 1267. 57i. However we will use that z is a zero in the following 
calculation. 

Applying the functional equation cos(^)r(s)C(s) = C(l — s)tt s 2 s ~ 1 and 
the reality of ( on R, 

(1 - z)((z - 1/2) _ (1 - z)(z - l/2)T(z - l/2)C(s - 1/2) 
C(z + 1/2) r(z + l/2)C'(z + l/2) 

(1 - z)(l/2 - z) (1 - sin(Trz)) C(3/2 - z) r 
2^ cos(ttz) C'(l/2 - z) ~ 

f.lzn 2/ 2 + 1/2 , / C(3/2 + z) y 

(5 - 14) —2^ { - l) Lav2 + z)) ' 

where z = iy, y S> 0. Thus we need to check that 

P (y) =^(C(^ + iy)(C(^ + iy)T) = »(<C(§ + %)(C^ + <»))*) = 

(5-15) »(C(| + n,)3(C*4 + - 9f(C(| + iy)»(C„(^ + ij/)) 

is mainly positive for ^(1/2 + iy) = <?£(l/2 + iy)/dy. The first term of the last 
difference dominates (we will skip the explanation). Let us show that it has a 
very good reson to be positive on the zeros z = iy. 

In the first place, 5ft(£(3/2+iy)) is always positive ([E], 6.6, pg.129). Then, 
following [E], 6.5, 

(5.16) C(\ + W) = e'^ y) Z(y), 6{y) = 9 log(r(i + i V -)) - | log(vr) 

for the classical real- valued function Z(y) of y £ R+. For instance, the zeros 
of 3 : C(l/2 + iy) = —Z{y) sin(0(y)) are either the zeros of Q(k + 1/2) (coming 
from Z) or the zeros of sin(0(y)), the so-called Gram points. For small y, they 
appear alternately (pg.125, Gram's law), i.e. indeed 9(^(1/2 + iy)) > on 
the imaginary (+0.5) zeros of the zeta. Even when Gram's law fails the rule 
due to Rosser et al. states that excluding the zeros from the so-called Gram 
blocks the sign is plus (8.4, pg.180). It holds for at least 13.4 millions of zeros 
of C, many more than we need here. 

However the above reduction and the last argument are not sufficient 
to explain why the first "negative" zero for the initial r/ + (z\2) (see (5.11)) is 
z = 1977.27141 The sign of ^s(C y (l/2 + iy)) becomes minus for the first time 
at y = 282.4651. The next one is y = 295.5733. Not very impressive. 

There is something remarkable about z = 1977.2714i. Exactly at this 
point, the computer program looking for the zeros of the zeta has to become 
more sophisticated. The reason is that this zero is very close to the previous 
one z = 1977.17391 If there is a correlation between the zeros with the linear 
^-approximations to the left and such pairs of zeros then it could be of a certain 
importance. The existence of nearly coincident zeros with low extrema of Z 
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between them "almost" contradicts the Riemann hypothesis ([E], 8.4, pg.178) 
and, indeed, "must give pause to even the most convinced believer" . 

6. Symmetrization 

In conclusion, let us try to establish connections with Riemann's esti- 
mate of the number of zeros in the critical strip and the Riemann hypothesis. 
We will examine the (anti)symmetrizations of the g-zeta functions considered 
above with respect to k <-> —k. It is the simplest way to ensure the functional 
equation. Let us begin with 3™(k\d) from (3.4) and its meromorphic con- 
tinuation 3+(f (fc|<i) to negative from (3.10). Since they behave differently 
for positive and negative one may hope to avoid unwanted zeros upon the 
symmetrization. 

Let us renormalize 3+^ (cf. (3.13)) as follows: 

Mk) m TT (W +2 *)(W +1 ) 

_0 (1 - q j+k+1 / 2 )(l - qi +k+1 ) ' 

(6.1) = imT q (k\d), T: 9 q {k\d) = f 7 (fc)3S(fc|d), 
lim a^+Md) = %+{k\d) = 7 k^+(k\d), m > 0, 
lim a k sT g 9 (k\d) = 3+(k\d), $lk < 0, for 

$+{k\d) = d- 1 ' 2 - k {l-2 1 ' 2 - k )T{k + l/2)C{k + l/2). 
Cf. (3.12), (3.13), and (3.23). For m > 0, we set 

(6.2) T+™(k\d) = f a k ^ g (k\d) + a- k 3 n + q 9 (-k\d), 

lim Z S r(k\d) = %*T(k\d) = f k^(3+(k\d) - 3+{-k\d)). 

a^co y 

The function 3+ g can be naturally extended to a symmetric (even) mero- 
morphic function defined for all k £ C: 3+ q m (—k\d) = 3+J"(&|<i). It has zeros 
at k € 2-iriaZ. Its poles modulo 2naiZ form (see (3.8)) the set 

(6.3) V = (-A) U A U {1/2 + Z}, A = {-2& - Z+} U {-2tf - Z+}. 

They are all simple for generic a. 

Given < n < 1 such that {k + iR + Z} n V = 0, 

lim a r ^ K 3+ m ( K ~ r + iu ) = a"V+> 

(6.4) where V+ = hm ZTSk) > 0, r e N. 

Here we employ Theorem 3.4 and (3.14) and use that (7/? _1 )(k — r + i/i) (see 
(3.13)) tends to zero as N 9 r — > — oo. 
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Because of the appearance of the a-factors the function 3+ 9 (k) is not 
27rai-per iodic. However we may assume that a 2nai = 1, i.e. 

cwlog(a) = M for M € Z. 

Then it becomes 2Ti-periodic for T d = ira and we can readily calculate the 
number of its zeros minus the number of poles modulo 2Ti in the limit of big 

— S'UTfL 

r. Moreover, because of the reality of 3 +(? (k) on the real axis, if z is a zero 
(a pole) so is 2Ti — z. Hence we can find this number for the half-period too. 

Corollary 6.1. Let k be as above, O r the number of zeros o/3+J™(&) 
in the strip 

S r d = {n-r <%lk <r - k, 0<Qk< T}, 

P r the number of poles. We calculate zeros(poles) with multiplicities and with 
the coefficient 1/2 if they belong to the boundary. Then 

(6.5) Noo d = O r -P r = -(log -) for N3r» 0. 

7T 7T 



Formula (6.5) can be considered as a certain counterpart of the celebrated 
formula (see [E, 6.7]) 

(6.6) N(T) = £(l og ^-l) + 0(logT) 

for the number of the zeros z of the zeta with < 9z < T (in the critical strip). 
The total number of imaginary zeros of 3+ (k\d) such that < 9z < T results 
from this formula modulo the Riemann hypothesis: 

(6.7) N(T\d) = -(logT - 1 - log(2d)) + O(logT). 

7T 

All zeros of £(k + 1/2) are those of 3+ m . However the latter may have extra 
zeros. According to our analysis their number is finite for d < 1 and is not 
bigger than (T/tt) log 2 for d > 1 For big they may appear only near the 
zeros of (1 — 2 1 / 2=tfe ). Anyway the non-imaginary zeros do not change (6.7) too 
much. 

Concerning 3 +g (k\d), first, it must have exactly one zero approaching 

— sy'tYi 

each pole from the set {±A; = l/2 + 2Z_|_} as q —> 1. Indeed, the limit 3+ {k\d) 
has neither poles nor zeros there. It is likely that the same holds for the 
remaining poles at {±k = —1/2 + 2Z + }. At least it is true for sufficiently big 
1 3F£Ar | because of the stable Z-periodicity of 2>+ g {k\d). Numerical experiments 
confirm this: the "symmetric" (/-zeros practically coincide with those from 
(3.20) as a = 1 = d. 
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Second, there is a tendency for the poles from A to be moving accompanied 
by the zeros. We cannot prove that this always happens. Anyway the number 
of (/-poles equals the number of g-zeros inside a closed curve in the form C a — 
y/aCi for a given C\ if 

\a k f™(k\d)\ > \a- k T+ 9 q {k\d)\ for k G C a , 

assuming that JftCi > and 3+^(&|c2) has no zeros inside C a . Generally speaking, 
it holds. However it is not clear at the moment whether this can be used for 
quantitative estimates. 

Anyway if it is true (or almost true) that the poles of 3+ g (k\d) from A are 
compensated by the zeros in the limit a — ► oo, then we come to the following 
conjecture. 

Conjecture 6.2. Assuming that d > 7r/(2e), the q- deformation z q of 
any imaginary zero z = iy(y > 0) of 3+ m (k\d), namely a zero of3+™(k\d) 
approaching z in the limit q — > 1, exists and remains in the horizontal strip 
< Qz q < y + C for all < q < 1 and a proper constant C > 0. The zeros of 
3+ q m {k\d) in the box 

/ 7TQ / 7T(2 

(6.8) K a (c) = {-2W — + c < m < 2J— - c, < %k < -no} 

y Ziili y Ziili 

are q- deformations of the classical ones for a constant c > 0. 

The first part of the conjecture for d = 7r/(2e) results in the estimate 

N(T) < £(log£-l) + 0(logT), 

speculating that the number of poles of 3+™(fc|<i) for big a with < < t 
which have more than one neighboring zero is O(logt) at most. As a — > 
(i.e. q = exp(— 1/a) — > 0), the g-deformations must approach the real line. So 
eventually they go down for small a. The box K a (0) in the second part does 
not contain the poles of Thus we expect that once a q-zevo gets distant 

enough from the g-poles it always approaches a classical zero in the limit. 

Variants. Now let us modify the imaginary g-zeta to avoid the gamma- 
factors. For Ifck > 0, we set 

( +q (k\d) = f Aa(k)3T q (k\d) - Aa(-k)3 n + e q 9 (-k\d), 

oo 

(6.9) where a(k) d = f J|(l - q k+1 ^)/(l - qi +1 ), 

j=0 

lim( +q (k\d) = ( + (k\d) = f (+(k\d)-(+(-k\d), 
C+(k\d) = (d/4)- 1 /2-fe(i_2 1 /2-fe )C(A: + 1/2) . 
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The function ( +q can be extended to an odd meromorphic 27rai-periodic func- 
tion with the poles at {irai + Z/2}. Since {a(5~ 1 )(k — r + v%) —>■ as r — > (see 
(6-4)), 

(6.10) lim ( +q (k\d)(K - r + iu) = Atp + , r G N. 

r^oo 

Hence the (stable) numbers of zeros of ( +q (k\d) and poles on the half-period 
< $sk < T = na coincide. So there is no relation to the Riemann estimate 
in the considered case, which can be explained as follows. First, the classical 
function £ + (k\d) has infinitely many real zeros approaching the points from 
±{1/2 + 2Z + } for big \$tk\. It can also have zeros which are neither real nor 
imaginary (but finitely many). Second, the poles of ( + (k\d) on the line = ira 
may "loose" the neighboring zeros as a — > oo. In line with Conjecture 6.2, such 
lost g-zeros could be a significant source of the classical ones. Hopefully it may 
be seen numerically. 

It is instructional to get a g-deformation of the modified zeta 

satisfying the plane functional equation ((k) = ((— k). As always, we need to 
switch to its plus- variant to ensure the convergence for Wt < 1/2: 

( + (k\d) = P (k\d)C(k), P {k\d) = 

(6.11) d-V 2 - k (l - 2 l ' 2 - k ) - d-Va+^i _ 2 1 /2+fc). 

All zeros of the function p(k\d) are imaginary for d > 1/2. The number of 
imaginary zeros in the interval < ^sk < T equals \{T/ir) log(2d)| up to —1 < 
e < 1 for any d > 0. Let us assume that d > 1/2. 
We set 

(6.12) a(fc) 7(*0 II x _ g fc/2 + i/4 + j ' (^K +q (k\d) = 

(^fl^a{k)^ q (k\d) - (^)- k ^M-k)37 q 9 (-k\d)- 

Then lim a ^ 00 C+g(fe|^) = (+(k\d). 

The function {a(3~ l )(n — r + vi)) tends to zero as r — > oo too, so we 
can calculate the stable (see above) difference between the total number 
of zeros and poles on the "half-period" for ( +q (k\d). We may assume that 
(a/(16vr)) 7rai / 2 = 1, i.e. 

(7rai)log(o/(167r)) = M for M € Z, 

to provide the 27rai-periodicity of ( +q . Then = M. Let us express it in 
terms of T = 7ra : 

(6.13) iVoo(T) = M = ^(log^-logvr). 
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Here the leading term coincides with that from (6.6) again. More exactly, the 
number of zeros of ( + (k\d) such that < Ssz < T equals 

(6.14) N(T\d) = ^(log£-l + 21og(2d)) + 0(logT). 

The zeros of ( +q either a) approach its real poles at {±k = 1/2 + 2Z + }, 
or b) go together with the corresponding poles from A (see (6.3)) upwards, or 
c) tend to the zeros of zeta in the critical strip. We may conjecture that a q- 
deformation z q of a zero z = iy(y > 0) of C(k + 1/2) remains in the horizontal 
strip < Qz q < y + Cdy/logy for all < q < 1 and a constant Cd > 
depending on d. One can also speculate that Cd ~ log(32-7r<i 2 /e) comparing 
the formulas for N(T\d) and N^T) in more detail. 

The above estimates can be extended to the case of the standard zeta. 
The factors (1 — 2 1 / 2 ~ fc ) will disappear from the formulas. The calculation 
of the stable number of zeros minus poles remains the same. However the 
minus-counterparts 3^"\ C?> an d C? °f the functions considered above diverge 
in the critical strip. So the relations to the classical zeros have to be discussed 
under the constraint log(a) = o(\T(k) 2 tan(7rfe)|~ 1 / 2 ). The conjecture and other 
considerations can be also extended to the (/-deformations of the Dirichlet L- 
function (see (5.4)), where the convergence holds for all k except for the poles. 

Towards Riemann g-hypothesis. The conjecture states nothing about 
the pure imaginary zeros of Z^+T (k\d). It is more convenient to discuss this 
problem for the sharp zeta-function. For the sake of concreteness, we assume 
that d = 1 and will (/-deform the ( from (6.9) only. A similar behavior is ex- 
pected for the sharp-variants of the remaining two symmetric q-zeta functions 
considered above, their minus-counterparts (for a which are not too big), and 
the sharp Dirichlet L- functions. 

We set 

(6.15) C[ q (k) d = f xiky^Uk) - xi-kr^U-k), 

oo 

JJ(1 _ _ + g j/2-l/4+fc/2 )(1 + q j/2-l/4-k/2^ 

3=1 

hmC+#) = C+(k) = d ^ C+(*0-C+(-*0- 

2 

Thus we just divided the sharp-integral for the kernel (1 + q~ x ) by that for 

2 

the Gaussian q x . This very choice of the normalization is not too important. 
However in this form the definition can be naturally extended to the multi- 
dimensional case. Thanks to Theorem 2.3 and Theorem 3.1, the function 
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C\ q (k) is regular in the horizontal strip 

Kl(c) d = f {-V2^ + c<%k< V2^ - c} 

for c > 0. 

Conjecture 6.3. Let us fix q = exp(— 1/a) for a > 0. Then all zeros 
of C+{k) inside if 1(0) have unique q- deformations, i.e. the zeros of C+(k) 
convergent to the corresponding classical ones in the limit, which belong to 
K+{—c) for a constant c > 0. Vice versa, all zeros of(^_(k) inside K\(c) are 
such deformations for proper c > 0. The q- deformations of the imaginary zeros 
°f C+(k) in K |(0) are imaginary. 

The sharp-counterpart of the Riemann hypothisis is the claim that all 
zeros z q of C+ g (fe) in the critical strip {\$lk\ < 1/2} such that < Qz q < 
\f 7 2/na — c for a certain constant c > are imaginary. Confirmations of the 
conjecture are entirely numerical. No theoretical approach is known at the 
moment. 
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